Comment on "Critique of the foundations of time-dependent density-functional 

theory" [Phys. Rev A. 75, 022513 (2007)]. 
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A recent paper (Phys. Rev A. 75, 022513 (2007)) challenges exact time-dependent density func- 
tional theory (TDDFT) on several grounds. We explain why these criticisms are either irrelevant 
or incorrect, and that TDDFT is both formally exact and predictive. 

PACS numbers: 31.15.Ew,71.15.Qe 



Time-dependent density functional theory (TDDFT) 
has a rigorous foundation P, i, i, H Q. Schirmer 
and Dreuw (SD) appear to criticise TDDFT on several 
grounds almost all of which are ultimately conceded 
by SD themselves. For example, SD state that "...an er- 
ror is introduced in both the TD and static KS linear 
response theory if the perturbing (external) potential is 
given by a nonlocal operator" , but this is misleading since 
such potentials do not exist in DFT. Another example is 
their long discussion of the problems of the RG action. 
Such problems were first raised a decade ago 0,01 j and 
resolved shortly after using the Keldysh action Q. The 
problems with the action in Ref. [l| have been thoroughly 
investigated in several works 0,1, Sii- Finally, SD 
admit these points, but by not addressing the current 
literature, they are tilting at windmills p^ . 

The sole unresolved criticism of SD is their claim that 
the Kohn-Sham (KS) equations of TDDFT, even if able 
to reproduce the density-evolution of the true interact- 
ing system, n(r,t), cannot predict that evolution. The 
KS potential, Vs (r, t) , is defined as the unique one-body 
potential in which non-interacting electrons evolve with 
density n(r, t). SD claim that their "radical KS" scheme 
shows that Vs (r, t) functionally depends on the future 
density, thereby making direct propagation impossible. 
They further argue that, with neither a variational prin- 
ciple nor some proof of convergence of the TD KS equa- 
tions for such a potential, TDDFT must be "unfounded" . 
The rest of our Comment addresses this issue: we show 
where the error lies in SD's arguments and explicitly 
demonstrate that the KS-TDDFT procedure is indeed 
predictive. 

The flaw in SD's arguments arises from not care- 
fully distinguishing between Vs{r,t) and the exchange- 
correlation (XC) potential. If v{r,t) is the time- 
dependent external potential, i.e., the one-body poten- 
tial applied to the interacting system (which is always 
known and given), and VH{r,t) = J n{r' ,t)/\r — r'\d^r' is 



the Hartree potential, then 
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is known to be a functional of the initial states (both 
interacting, 5*0, and KS, $0) and n{r,t), written 
Wxci"-, ^'o, ^o](r, f). SD argue (pg 11) that, since Vs{r,t) 
depends on dfn{r,t) (the second time-derivative, which 
SD denote as "the future"), the possibility of a stringent 
time propagation is "thwarted" : if the potential func- 
tional depends on the "future" of the density, then how 
could it also possibly predict it? In particular, at the 
initial time t = Q, n{r, 0), ^"0, and <I>o are obviously insuf- 
ficient to determine dfn{r,0), and therefore insufficient 
to determine Ws(r,0). 

What SD miss is that in any TDKS time propagation, 
w(r, t) is known and prescribed by the physical problem 
at hand, and it is only Wxc(r, t) for which a density- initial- 
state functional is needed. Herein lies the resolution of 
the propagation paradox, as the functional dependence of 
Vs{r, t) generally differs from that of Wxc(r, t), in contrast 
to what is claimed by SD without justification: "A sim- 
ilar temporal nonlocality must be expected for the XC 
potential...". In particular, Wxc["-, *o, 'i'oKr, 0) depends 
purely on the initial states, as we explicitly show below, 
and SD's expected dependence on d^n{r, i = 0) vanishes. 
At other times, unlike v{Y,t) and Vs{Y,t), Wxc(r,^) does 
not depend on dfn(v, t): there is a dependence on dtn{rt) 
but this can be extracted from $(t), through the conti- 
nuity equation, and so is available during propagation. 
(There is also dependence on 4'o and (f>o and the earlier 
history of the density). 

Before demonstrating in general the density- and 
initial-state-dependence of Wxc(r, t) ("potential function- 
al" in SD's notation) and why direct propagation works, 
we give a simple counterexample to SD. Propagate the 
time-dependent Schrodinger equation for one electron in 
some potential from some initial state. Applying SD's 
"radical KS" scheme and logic, this is impossible, be- 
cause that potential functionally depends on the future 
density (eg. Eq.75 of SD). This absurd conclusion is in- 
correct, because you are given v(r, t). You never extract 
it from the density. Instead, what is needed for standard 
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KS TDDFT is a dcnsity-initial-state functional only for 
Wxc(r,i)- For one electron, t)xc(r,i) — ~v^{v,€) (choos- 
ing the KS initial state as the true initial state, the only 
sensible choice [ll|), i.e. a functional of the instantaneous 
density alone. There is no dependence on dfn{r,t), con- 
trary to the "expectations" of SD. More generally, be- 
cause only the density-dependence of Wxc(r, i) is needed 
in any TDDFT KS calculation, their result for Ws(r, t) is 
irrelevant to the question of propagation. 

First, to clarify the question of "future" dependence, 
we note that a dependence on dfn{r,t) does not violate 
causality at any finite t. For any < t <T, v[n; '^o]{r, t) 
can depend only on the density in the interval [0, T] and 
not on later times because, by the RG theorem, two po- 
tentials that are the same in [0, T] but differ ior t > T 
must give different densities at times t > T, i.e., v is 
a causal functional of the density. The same holds for 
Vs[n\^o]{Y,t). From Eq. (1), Uxc[n; *o, *o](i) is then 
also a causal functional of the density and initial states 
for any finite t, and dfn{r, t) may be evaluated to the left 
of time t, (e.g. using n(<), n(t — Ai), n(t — 2At)). But this 
argument cannot be applied at t = 0, and so the start of 
a propagation appears problematic. 

To see why this is in fact not a problem, we will use 
results from Ref. [5], and the notation that a super- 
script (fc) indicates the fc-th time-derivative at < = 0, eg. 
rS^\r) = dfn{rt)\t=o- Begin by noting that both n(r,0) 
and n'^^^(r) are determined by 5*0, the initial wavefunc- 
tion of the interacting problem, because of continuity: 
dtn{r,t) — —V • j(r, t). This in turn restricts allowed 
choices of $0 to only those that recover these values. 
The form of v{r,t) influences only second and higher 
time-derivatives of n(r, t). The TD Schrodinger equation 
implies [H| 

d^n{r, t)=V- [n(r, t)Vv{r, t)] + q{r, t), (2) 

where 

q{r,t)^{^it)\ T{r)+w{r) \^it)), (3) 

and 

= ^^a.^fe (9,^t(i.)a,^(r)-f 5fe^t(i.)a,v3(r) 

i,k 

w{^) = ydVVt(r)v}t(r')9fc^--l-^^(r')^(r). 

An analogous equation applies to the KS system, with 
Vs{rt) and qs{r,t) = {^{t)\ f(r) |$(i)}. Requiring that 
the density evolutions be the same yields 

V • [n(r, t)VvHxc(r, t)] = q{r, t) - qs{r, t) (4) 

where WHxc(r, t) = Vii{v, i) -|-Wxc(r, t). This equation is of 
Sturm-Liouville form: such equations have a unique solu- 
tion for UHxc(r, t) if n(r, t) and g(r, t) — qs{r, t) are given. 



together with the boundary condition that Whxc(^ ^ 
00, t) 0. We shall assume this boundary condition 
for all that follows; any choice of a TD constant for the 
asymptotic potential does not affect the density [5]. Thus 
a KS potential can always be found for any density for 
the interacting system, provided the initial conditions are 
met. 

We now use Eq. @ to show that at each time step, the 
functional input to Wxc(r, consists of the ^0, 'i'o: and 
the density evolved through previous times. At f = 0, 

V • [n(r, 0) V«Hxc(r, 0)] = (/(r, 0) - (Zs(r, 0) . (5) 

Since q{T, 0) and qs{Y, 0) are determined by Eq. ([3]), which 
is known entirely from the initial wavefunctions, the solu- 
tion of this equation determines Wxc(r, 0). (The Hartree 
potential UH(r,0) is determined as usual directly from 
the instantaneous density). Thus, as mentioned ear- 
lier, Uxc(r,0) depends on the initial states alone: the 
dependence on the second-time-derivative of the density 
in Eq. ^ and its KS counterpart cancel once their dif- 
ference is taken in Whxc["-, ^Oi *i'o] (rO). Together with 
w(r, 0), this evolves $0 forward one time step: 

(^-iv2 + t-Hxc(r,0) + t;(r,0)^ (/),(r, 0) = ^at0,(r, 0) 

(6) 

This yields the orbitals at the first time-step, (/)i(r. At) — 
0i(r, 0) -I- 9t(/)i(r, 0)Ai, from which one can obtain 
the density evolved to the first time-step, n(r. At) — 
|0i(r, At)|^, as well as the current at the first time- 
step, j(r,At) = 3^^0*(r,At)V0,(r,At), and, through 
the equation of continuity, also the first-time derivative 
of the density rj(r. At) = -V • j(r. At). 

Next, consider t = At, for which we need to find 
(?t'f^Hxc(r, 0). Take a time-derivative of Eq. (jH), evalu- 
ating the terms at t = 0: 

V • [n(r, 0)V9tWHxc(r, 0)] = -V • [dtn{r, 0)Vi;Hxc(r, 0)] 

+at(g(r,0)-gs(r,0)) (7) 

where 

9t(g(r,0)-gs(r,0)) - *(<i>o|[f (r), 7?s(0)]|$o) - 

*(*o|[g(r),i?(0)]|vl/o) (8) 

Equation ([7]) is again of Sturm-Liouville form, with a 
unique solution for 9tt;Hxc(r, 0): all the other quantities 
in the equation are known either from the initial states 
or the results of previous time-propagation. Eq. ([7]) is 
not an explicit functional of the initial states and density 
alone, due to the appearance of w(0) in the commutators 
in Eq. ([5]) . But evaluating Eq. ([J) at t = gives the func- 
tional dependence of u(0) on 5'o,n(r,0), and 9j^n(r, 0). 
Since ^^^^(O) = (9tn(At) - 9tn(0))/At (limit At ^ 
understood), the appearance of the second-derivative at 
t = implies only a dependence on the first derivative at 
t = At, directly available from <i>(At) via the equation of 
continuity. So 

ws(r, At) = t;(r, At) + WHxc(r, 0) -f 9tWHxc(r, 0)At (9) 
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is determined, and predicts the time-evolution of the den- 
sity at t — 2A. 

Each subsequent time-derivative of Eq. ^ produces 
one higher time-derivative of the XC potential, such that 

^'Hxc(^) is determined solely by the initial states and 
n(')(r) with / < (fc-|- 1), all available from propagation to 
the fcth time-step: 



V-[n(r,0)V^;Sc(r)]-Q(''^)(r) 



(10) 



for practical applications [20|. The equations for these 
approximate potentials also show the fundamental prop- 
erty that we demonstrated above: the XC potentials at 
a given time are completely determined by the density 
evolved up to and including the present time and the 
initial states. 

Lastly, we give an example to show explicitly the error 
in SD for an interacting case. Consider two electrons in 
one dimension, and we need look only at t=0. In one- 
dimension, Eq. ([5]) reduces to 



QW(i.) = g(fe)(r)-q«(r)- 



fc-i 

E 

/=0 



V-[n('=-')(r)VwSc(r)] .Ws(a;, 0) = Wcxt(a;, 0) + 



dx' 



n{x',0) 



dx"{q{x,Oyqs{x,0)) 



(11) 

The q^'^^ (r) and qi'^'' (r) involve multiple commutators of 
the operator q(r) with the true and KS Hamiltonians, 
respectively, and their time-derivatives, sandwiched be- 
tween the mitiaZ states ^"0, *&0j respectively. For example, 
for fc = 2, 



But inversion of the time-dependent KS equation [l^[lc 
as in SD, yields here 



Vsix, 0) 



g(2)(r)-gf)(r) = 
-(vI/o|[(z(r), i?(0)], iJ(0)]|*o)+ (<fo|[r(r), 7?s(0)], iJs(0)]|<i>o 
-i{^o\[qir),dH/dm]\^o)+ i($o|[f(r),dFsM(0)]|$o) 

The v{0) and dv{Q)/dt appearing in the commutators are 
causal density and initial-state functionals via Eq. ([2]), so 
Eqs. (fTU]) and pT|) yield the XC potential as a causal im- 
plicit functional of purely the density and initial states. 
We reiterate that a dependence on dtn{r,t) presents no 
difficulty during propagation. As one propagates, one 
can evaluate this instantaneously, by simply computing 
the divergence of the current-density of ^{t). Thus all 
quantities needed are available from the past propaga- 
tion. 

So we have shown that the propagation can be done 
in a predictive manner, and that the expressions for 
Vxcin, Vfo, 'i'oKr, t) that we give are causal functionals of 
the density and initial-state: explicit at t=0 (Eq. and 
implicit at later times (via *o], Eqs. l4l7l8fTT|) . This 
construction holds for any time-dependent potentials and 
densities that are equal to their Taylor expansions for 
i > for a finite period of time (as assumed for v{r, t) 
in the Runge-Gross one-to-one mapping proof). Full an- 
alyticity is not required: in particular, as nothing is as- 
sumed for times earlier than t = 0, this procedure applies 
to sudden switch-on potentials. 

The XC potentials used in practical applications may 
be viewed as approximations to this formally exact con- 
struction of the potential functional. Although most of 
the applications to date have utilized adiabatic approxi- 
mations, depending on the instantaneous density alone, 
memory-dependence is a well-recognized feature of time- 
dependent functionals (see eg. Refs. [H, [H, [H, [H, 
|l9[) and memory-dependent functionals have been suc- 
cessfully applied to real calculations. There are also ex- 
plicit systematic methods based on many-body pertur- 
bation theory to construct approximate XC potentials 



dxTl 

dtj 
n 



2 V" 
2^2 



in 



dx' 



(13) 



)where j = j{x,Q) is the initial current-density, deter- 
mined from the initial wavefunction. (We drop the spa- 
tial and t = indices on the right, for ease of reading). 
The alleged dependence on the future arises through the 
term dtj on the right: this may be equivalently written 
in terms of dfn{x,0), and is the only term not directly 
obtainable from the initial states. 

However, this future dependence disappears as soon as 
we relate the KS system to the interacting system via 
the Heisenberg equation of motion for the current of the 
interacting system: 



dtj = -ndxVc 



Tx. 



Wx 



(14) 



where d^T^x = (*o|t|*o) and d^W^ = (^'o|i&|*o)- Sub- 
stituting this for dtj in the right-hand-side of Eq. 



^dx' = 



dx' 



n{x',0) 

= Vcxt{x,0) + 



d^n{x",0)dx" 
Txx + W, 



dx' (15) 



i.e. the "future" dependence is in the external potential 
and other terms at i = 0! That is, once the connection 
with the interacting system is made, the apparent depen- 
dence on the future evaporates as d1n{x^ 0) is determined 
by initial-state information and by the external potential 
that the interacting system is subjected to. The depen- 
dence is explicit in this two-electron example, but in the 
general A^-electron case, the construction of Ref. Q im- 
plies that this is always true. 

Inversion of the KS equations alone yields information 
only about the KS potential as a functional of the den- 
sity, but tells nothing about the XC potential. How could 
it, since it contains no information about any interacting 
system? This is most easily seen in the ground-state 
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problem. One can trivially invert the single KS orbital 
equation for any two-electron density and get its KS po- 
tential, but this tells you nothing about XC unless you 
know the corresponding external potential to subtract 
from it. And there's no way to find that, without invert- 
ing the interacting Schrodinger equation, thereby making 
the functional dependence as implicit as in the original 
definition fl^,^. 

Before concluding, we briefly review SD's discussion of 
the lack of a proof of numerical convergence of KS propa- 
gation when non-adiabatic functionals are used (the "tra- 
jectory mode" in SD's notation). SD only claim that such 
a proof is needed after they incorrectly deduce that the 
"potential functional" propagation mode fails. This is 
therefore a straw man argument. There are no proofs of 
the existence of solutions and their convergence for, e.g., 
the Navier-Stokes equations, but should they be "aban- 
doned?" (In fact, the KS propagation has been explicitly 
converged in a recent calculation with a non-adiabatic 
potential 

To summarize, (i) SD's rejection of the original RG 
formulation of TDDFT originates in an elementary log- 
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